Anholonomic spin manipulation in drift transport in semiconductors 
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We find that the electronic spin rotation induced by drift transport around a closed path in a 
wide variety of nonmagnetic semiconductors at zero magnetic field depends solely on the physical 
path taken. Physical paths that produce any possible spin rotation due to transport around a closed 
path are constructed for electrons experiencing strain or electric fields in (001), (110), or (lll)-grown 
zincblende semiconductor quantum wells. Spin decoherence due to travel along the path is negligible 
compared to the background spin decoherence rate. The small size of the designed paths (< 100 nm 
scale in GaAs) may lead to applications in nanoscale spintronic circuits. 
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Anholonomy, in which a vector transported around a 
closed curve in parameter space does not return to its 
initial orientation, is a general property of quantum me- 
chanical systems, and has analogs in other fields includ- 
ing optics and classical mechanics. [1] Within quantum 
mechanics the geometrical anholonomic phase, which de- 
pends only on the path taken and not on the time taken 
to traverse the path, is the Berry phase. [2] Recently the 
spin precession during drift transport of polarized pack- 
ets of electron spin in semiconductors has been shown to 
be independent of transit time when no magnetic field 
is applied. [3-6] The source of this behavior is the special 
nature of the spin-orbit Hamiltonian in these systems, 
where the lack of inversion symmetry combined with 
strain[4-6] or confinement in a quantum well[3] produces 
a spin splitting that depends linearly on the electron crys- 
tal momentum. Spin precession that is independent of 
transit time, however, is not sufhcient to demonstrate 
anholonomic transport and a geometrical phase, for re- 
versal of the path taken will unwind the spin precession 
such that the electron spin will return to its original ori- 
entation. Anholonomic manipulation of spin in quantum 
dots has been described [7-9], however this effect requires 
coherent manipulation of the confined electron's wave- 
function and does not apply to spin manipulation of an 
orbitally-incoherent ensemble of spins undergoing drift 
transport. 

Here we demonstrate that the drift of electron-spin 
polarization in a semiconductor, under the influence of 
strain or quantum well confinement, corresponds to a 
general nonholonomic system, whereby closed paths can 
be constructed that will produce any desired spin pre- 
cession on an initial spin, depending only on the path 
and not on the transit time. Ordinarily in an anholo- 
nomic system the Berry phase must be separated from 
the dynamical phase, which does depend on the transit 
time; in these nonmagnetic semiconductor systems the 
dynamical phase vanishes, leaving only the Berry phase 
(an SU(2) rotation of the spin) behind. Corrections to 
the linear dependence of the spin splitting on the electron 
crystal momentum produce dynamical phases, and also 
spin decoherence. We identify realistic situations where 
these corrections are small. Typical scales of the paths 



in 10 nm quantum wells are 50 — 2500 nm, leading to 
the potential for sub-100 nm scale spin manipulation el- 
ements. Anholonomic spin phenomena will permit the 
preparation of arbitrarily-oriented spin packets for prob- 
ing fundamental spin dynamics as well as new methods of 
electrically controlling spin resonance. Current interest 
in spin manipulation for information processing[10, 11] 
may motivate the use of this effect for spintronic devices. 
Spin anholonomy is possible with any of the typi- 
cal spin-orbit fields generated in semiconductor quan- 
tum wells, including the so-called Rashba, Dresselhaus 
and strain-induced fields, and for (001), (110), and (111) 
quantum wells. Each of these spin-orbit interactions can 
be described by a precession field r2(k) = — 0(— k) which 
depends linearly on the crystal momentum k, and the 
Hamiltonian 

Hso = S-n{k)='^a-n{k)='^a-E-k (1) 

for the electrons of the lowest-energy conduction band, 
where 2 is a second-rank tensor (the left-side dot prod- 
uct in Eq. (1) involving one set of indices and the right- 
side dot product the other) . The effect of a Hamiltonian 
such as Eq. (1) on spin precession in transport has been 
observed[4, 5, 12, 13] and utilized to initialize and manip- 
ulate spins through control of the transport direction[14]. 
The correspondence between precession lengths in trans- 
port and forms of the Hamiltonian has been extensively 
investigated [15-17]. 

We describe the motion of a population of carri- 
ers with uniform total density, but a spatially and 
temporally-dependent spin polarization using the semi- 
classical Boltzmann equation[18]: 
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The condition of uniform total density can be enforced 
by introducing an Ansatz for the spin-dependent density 
matrix F, such that F(x, k, i) = (I -t- m(x, t) • CT)/(k), 
where the trace of /(k) over k yields the uniform den- 
sity no, I is the identity matrix, and m(x, i) is the local 
(dimensionless) magnetization. We then make the relax- 
ation time approximation, introducing I[/(k) — /o(k)]/T 



for the spin-independent portion of the coUision term 
on the right hand side of Eq. (2) (where /o(k) is the 
equihbrium orbital occupation function and r is the 
orbital relaxation time) and (m(x, i) • a)f{'k)/Ti for 
the (isotropic) spin-dependent portion of the collision 
term. Anisotropic spin relaxation can be described with 
anisotropic relaxation times Ti for the different magne- 
tization directions. When this Ansatz is introduced, and 
the resulting equation is traced over spin, the ordinary 
spin-independent steady-state Boltzmann equation is ob- 
tained, (el^/h) ■ Vk/(k) = -[/(k) - /o(k)]/T. 

By taking the trace of Eq. (2) with different Pauli ma- 
trices, and using the definition of the (uniform) charge 
current jo = nge/iE = Trk(e/ik/7n*)/(k), we find 



9m(x, t) 
dt 



-Ai(E-Vx)m(x,i)- 



-^m(x,i)xS.E-"^(^ 

(3) 
We note two special cases of particular relevance to spin 
transport. If the spin polarization is uniform in space, 
then the Vxni(x, t) term on the right hand side vanishes 
and we find a direct relationship between the electric 
field and the temporally-uniform precession rate of the 
spatially-uniform spin polarization. For our second case, 
the steady-state solution, the left hand side vanishes and 
in the absence of spin relaxation we obtain 



(E • Vx)m(x, t) + -^m(x. f) x S 
n 



0, 



(4) 



where E = E/_E. Eq. (4) describes propagation that is 
independent of time and of the strength of the electric 
field (only depending on the direction of the field) . Fur- 
thermore, the gradient of m is perpendicular to m, so (as 
we have neglected spin relaxation), the magnetization's 
magnitude does not change with time. 

Eq. (4) is a central result of this work, for it shows that 
the evolution of a local magnetization of spin-polarized 
electrons can be described fully by the evolution along 
a path of the position of those electrons. If the electric 
field direction is changed with time, Eq. (4) will hold 
except during those transient periods during which the 
electric field direction is changing or the occupation func- 
tion /(k) is changing to its new steady-state value. For 
an abrupt change in E, the timescale of this change will 
be the orbital relaxation time t, so for Eq. (4) to dom- 
inate the spin transport, r must be much less than the 
length of time the electric field is constant as the spins 
are transported around a path. 

The precise form of S originates from a variety of phys- 
ical configurations of the semiconductor, including strain 
and quantum- well confinement. We use a notation capa- 
ble of describing all (fc-linear) precession fields in (001), 
(110), and (111) quantum wells, with their accompanying 
possible as-grown strain: 

hrt{]i.) = a{kyTi - k^y) + /3{kx:X - kyf) + r]kyi. (5) 



Here a is the Rashba coefficient (for (111) wells there 
are other contributions to a). In an unstrained (001) 
quantum well, j3 is the linear Dresselhaus coefficient and 
ry = 0, whereas for (110) growth 77 is the linear Dres- 
selhaus coefficient and /3 = 0, and for (111) growth 
r] — (3 = Q. Strain from growth on a lattice-mismatched 
subtrate modifies the values of a, /3, and 77. It is conve- 
nient to express Eq. (5) in a form that will lend itself to 
the path description developed below, 
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(6) 

where k? — a^ + P'^ , tanp = a//3, z is perpendicular to 
the plane of the quantum well, and k and 4> are the mag- 
nitude and azimuthal angle of the electron momentum 
(confined to the quantum well plane). S is then 
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For example, when a = 77 = in Eq. (5), Vtti^: ~ 
{m* P/h'^)mzy, ^my ~ {m* f3/fi^)mzX, and Vm^ — 
-{m*l3/h^)(m^y + my±). 

In general, motion along a straight path yields a spin 
rotation determined by the direction of the path, 0, and 
the distance travelled, r. To simplify the description we 
consider 7] — Q and find from Eq. (4) the electron spin 
rotation angle 



m Kr , ,1/0 

X= —p^il + am2(t)sm2py'^ , 



(8) 



and the (in-plane) azimuthal angle of the precession axis 



$ = - tan 



sm{p + <p) 
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(9) 



The same type of result holds for 77 7^ and /? = 0, 
although the precession axis lies in a plane oblique to 
the growth plane. The length for the spin to precess 2tt 
(precession length). 



^p{(t>) = ^^(l + sin2(/)sin2p)-i/^ 



(10) 



In both the pure Dresselhaus {p — 0) and the pure 
Rashba {p = 7r/2) cases the precession is isotropic, thus 
X = m*Kr/h? and Ao((/)) = ^.■K/2{4') = 2-Kh^ /{rrfn). For 
other p the precession length varies with the direction of 
travel (</)) in the quantum well plane. 

The spin rotation when traveling along a path can be 
constructed by sequentially solving Eq. (4), if the change 
in /(k) at the points between line segments of the path 
does not lead to additional spin precession or spin relax- 
ation. This condition can be met if the orbital relaxation 
time T is much shorter than the path transit time (which 
is also a condition of the relaxation time approximation) . 



We describe the motion around a path as a product of 
spin precession operators C/„ . . . C/2 C^i , where each 
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(11) 
that cor- 



arises from the displacement r with azimuth 
responds to one segment of the path. 

We now describe path constructions to generate an ar- 
bitrary spin rotation about any axis (the operator i?(0)), 
where the desired precession vector Jl is described by x-> 
$ and 9. We focus on closed paths, because it is pos- 
sible to transport a spin an arbitrary distance without 
precession by combining line segments made up of whole 
numbers of precession lengths. We find that a path with 
at least three non-collinear sides i?(ri) = U3U2U1 can 
produce a general rotation without net spatial displace- 
ment, although the solution is not unique. 

We begin with a rotation about an in-plane axis. This 
rotation can be made with a single line segment Ui , how- 
ever to achieve spin manipulation with zero net displace- 
ment it is necessary return to the origin by a different 
route, with each leg causing a 2tt rotation. For rotation 
axes in the quantum well plane (Q = 0), constraining 
U3 = U2 = I (by making them A in length) and minimiz- 
ing the total length will yield a unique solution, corre- 
sponding to the triangle shown in Fig. 1. The principal 
leg of the triangle is of length £ and azimuth 0, which are 
derived from x and $ by inverting Eqs. (8-9): 



h^X 



m* K\/l + sin 2p sin 20 
_^ ( cos (p — $) 



(12) 



VI + sin 2/9 sin 2$ 



The remaining vertex of the triangle is placed at an in- 
tersection of the ellipses of constant phase (with radius 
A((/)')) centered at each end of the principal leg. The sec- 
ond and third legs of the triangle return the spin packet 
to the origin having performed only 27r precessions — 
equivalent to identity operations — about arbitrary axes. 
Controlled rotation of the spin about an arbitrary 
(non-planar) axis is also possible. The precession field 
is always confined to a plane (the quantum well plane in 
(001) and (111) wells, but an oblique plane for (110) wells 
with a Rashba field), however successive non-collinear 
displacements result in non-commuting rotations of the 
spin, which can generate out-of-plane and arbitrary direc- 
tion rotations. As a special case, for coherent transport 
(not drift transport as considered here), a square loop in 
(001) quantum wells causes rotations proportional to the 
side length about the (001) axis[17, 19, 20]. An arbitrary 
rotation axis can, in principal, be implemented by a tri- 
angular path, but it is more intuitively implemented by 
the five sided path shown in Fig. 2. The principal seg- 
ment, of length £, is bracketed by two legs of length Ap/4; 
the first rotates the precession axis il into the plane of 





FIG. 1. (color online) Sketch of a path which will rotate the 
spin polarization of an electron about any in-plane axis il with 
azimuthal angle $. Unlabeled legs are nAp in length, where 
n is an integer, (f> is the azimuthal angle of the principal leg. 




FIG. 2. (color online) Sketch of a path which will rotate the 
spin polarization of an electron about any arbitrary axis f2 
with azimuthal angle "1> and elevation O. The projected axis 
of rotation il\\ with azimuth $ is also shown. The auxil- 
iary legs have length Ap((^)/4 and azimuth </> from Eq. (12). 
The principal leg has length I and azimuth 0' calculated from 
Eq. (12) using $' = "1> — O in place of $. Unlabeled legs are 
nAp in length. 



possible rotation axes and the second restores it to its 
original orientation. The remaining legs of the path then 
return the packet to the origin with identity operations 
on the spin. 

Use of such paths in devices offer great flexibility. For 
example, by adjusting the Rashba field with an applied 
electric field so that p = 7r/4 (corresponding to the spin 
helix condition [21]) any motion of a spin around a closed 
path yields no spin manipulation. Changing the Rashba 
field away from p = 7r/4 provides spin manipulation — 
permitting gated on-off control of the spin manipulation. 

We continue with specific descriptions of the quantities 
a, /3 and 77 for quantum wells and strained semiconduc- 
tors, and estimates of the size of these closed paths for 
several materials and quantum wells (tabulated in Table 
I). The area enclosed by the path will typically be near 
to or less than the square of the precession length, which 
is inversely proportional to the strength of the precession 
fields. There are two strain-induced precession fields, one 
dependent on the shear strain and the other on the biax- 



ial. They are determined by 

hilx = (~J3\f^y^xy ~ rtz^xz) + (-^Akxi^yy ^ &zz) (I'j) 

and the other components are obtained by cychc permu- 
tation. We now describe the spin precession fields for a 
quantum well strained by growth on a lattice-mismatched 
substrate. In both the (001) and (110) quantum wells a 
is the Rashba coefhcient. For (001) growth rj = Q and 



/3 = -7(fc|)+C4 



^^xxxx r ■^^xxu-y ^tt 



(14) 



where Cijki and a are the stiffness tensor and lattice con- 
stant for the well material. For (110) growth /3 = and 



zxxx I ^^xxyy 



r; = 7(fc2)/2+(C3-C4) """" ' """^^ — . (15) 

^xxxx \ Cxxyy ~r Cxyxy ^ 



For (111) growth f3 = tj — {) and 



-^Cxxyy lAQ, 



a = afl+— 7(fc2)+V3C3- ^^ 

V O ^xxxx ^^ ^^-xxyy y^ ^^xyxy ^ 

(16) 
where aj^ is the Rashba coefficient. 



w z 


GaAs InAs GaSb InSb ZnSe CdSe GaP InP 


5 nm (001) 


66 37 


1.1 


1.2 390 660 350 24 


10 nm (001) 


260 150 


43 


47 1600 2700 1400 970 


10 nm (110) 


530 290 


86 


93 3100 5300 2800 1900 


10 nm (111) 


230 130 


37 


40 1300 2300 1200 840 



TABLE I. The precession lengths (in nm) of a series of quan- 
tum wells of differing thickness and growth direction com- 
posed of various fff-V and If- VI semiconductors considering 
only the Dresselhaus field, w is the depth of the quantum 
well, and z is the growth orientation for it. 

Our concluding concern is the role of spin relax- 
ation on the spin transport phenomena described here. 
Spin relaxation in semiconductors occurs via differential 
precession of spins in the spin-orbit field (dcphasing), 
followed by spin-independent scattering that produces 
decoherence[10], and the relaxation rate is proportional 
to the electron mobility. In the drift regime the drift 
momentum of the Fermi sea is much smaller than the 
typical momentum of relevant states (Fermi momentum 
in the degenerate state, thermal momentum in the non- 
degenerate state), so the relaxation rate is not affected 
by transport of carriers (unless the transport leads to 
an elevated temperature, for which case the relaxation 
time at the elevated carrier temperature is the relevant 
quantity[22]). The transport time around the path does 
depend on the mobility, but also on the electric field, so 
it should be possible with sufficiently large electric fields 
to move spins around the path with negligible spin deco- 
herence. For example, the spin relaxation rate for a clean 
(/i '^ 5 X 10"^ cm^/Vs), 5 nm thick GaAs quantum well at 



room temperature is 20 ps[23], which in an electric field of 
10^ V/cm corresponds to a 10 /zm spin transport length, 
more than 50 times the path length required to perform 
a spin operation. The error rate in the spin polarization 
after traversing this path would thus be < 2%. 

We have proposed a set of structures that allow for 
arbitrary manipulation of the spin orientation of an elec- 
tron packet moving via drift transport, without magnetic 
fields or net displacements, in a variety of quantum well 
structures and geometries. The overall size of the paths 
the spin packets must traverse can be quite small, al- 
though it remains an open question whether this form 
of spin anholonomy will persist in the ballistic transport 
regime. 
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